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Abstract 

Search problems at various scales involve a searcher, be it a molecule before reaction or a foraging 
animal, which performs an intermittent motion. Here we analyze a generic model based on such type 
of intermittent motion, in which the searcher alternates phases of slow motion allowing detection, 
and phases of fast motion without detection. We present full and systematic results for different 
modeling hypotheses of the detection mechanism in space dimension 1, 2 and 3. Our study completes 
and extends the results of our recent letter [Loverdo et al. Nature Physics 4, 134 (2008)] and gives 
the necessary calculation details. In addition, a new modeling of the detection phase is presented. 
We show that the mean target detection time can be minimized as a function of the mean duration 
of each phase in dimension 1, 2 and 3. Importantly, this optimal strategy does not depend on the 
details of the modeling of the slow detection phase, which shows the robustness of our results. We 
believe that this systematic analysis can be used as a basis to study quantitatively various real 
search problems involving intermittent behaviors. 

PACS numbers: 
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I. INTRODUCTION 



Search problems, involving a "searcher" and a "target", pop up in a wide range of do- 
mains. They have been subject of intense work of modeling in situations as various as 
castaway rescu e fora g ing animals or proteins reacting on a specific 



DNA sequence 
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gies, the so-called intermittent strategies - combining "slow" phases enabling detection of 
the targets and "rapid" phases during which the searcher is unable to detect the targets - 
have been proved to be relevant at various scales. 

Indeed, at the macroscopic scale, numerous animal species have been reported to perform 
such kind of intermittent motion [20I . I21I while searching either for food, shelter or mate. 
In the case of an exploratory behavior, i.e. when the searcher has no previous knowledge or 
"mental map" of the location of targets, trajectories can be considered as random. Actu- 
ally, the observed search trajectories are often described as a sequence of ballistic segments 



interrupted by much slower phases. These slow, and sometimes even immobile 



2l| . phases 



are not always well characterized, but it seems clear that they are aimed at sensing the en- 



vironment and trying to detect the targets 



22| . On the other hand, during the fast moving 



phases, perception is generally degraded so that the detection is very unlikely. An example 
of such intermittent behavior is given by the C.elegans worm, which alternates between 
a fast and almost straight displacement ("roaming") and a much more sinuous and slower 
trajectory ("dwelling") [23]. During this last phase, the worm's head, bearing most of its 
sensory organs, moves and touches the surface nearby. 

Intermittent strategies are actually also relevant at the microscopic scale, as exemplified 
by reaction kinetics in biological cells Q, [24 1. As the cellular environment is intrinsically 
out of equilibrium, the transport of a given tracer particle which has to react with a target 
molecule cannot be described as mere thermal diffusion : if the tracer particle can indeed 
diffuse freely in the medium, it also intermittently binds and unbinds to motor proteins, 



which per: 
filaments 



"orm an active ballistic motion powered by ATP hydrolysis along cytoskeletal 
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for 



nt trajectories of reactive particles are observed 
example in the case of vesicles before they react with their target membrane proteins 
In that case, targets are not accessible during the ballistic phases when the vesicle is bound 
to motors, but only during the free diffusive phase. 



26|. 
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As illustrated in the previous examples the search time is often a limiting quantity whose 
optimization can be very beneficial for the system - be it an animal or a single cell. In the case 
of intermittent search strategies, the minimization of the search time can be qualitatively 
discussed : on the one hand, the fast but non-reactive phases can appear as a waste of 
time, since they do not give any chance of target detection. On the other hand, such fast 
phases can provide an efficient way to relocate and explore space. This puts forward the 
following questions : is it beneficial for the search to perform such fast but non reactive 
phases? Is it possible, by properly tuning the kinetic parameters of trajectories (such as 



the durations of each of the two phases) to minimize the search time? 



hese questions 
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have been addressed quantitatively on specific examples in 
where it was shown that intermittent search strategies can be optimized. In this paper, we 
perform a systematic analytical study of intermittent search strategies in dimensions 1,2 and 
3 and fully characterize the optimal regimes. This study completes our previous works, and 
in particular our recent letter j^l, by providing all calculations details and specifying the 
validity domains of our approach. It also presents new cases of potential relevance to model 
real search problems. Overall, this systematic approach allows us to identify robust features 
of intermittent search strategies. In particular, the slow phase that enables detection is 
often hard to characterize experimentally. Here we propose and study 3 distinct modelings 
for this phase, which allows us to assess to which extent our results are robust and model 
dependent. Our analysis covers in details intermittent search problems in dimension 1, 2 
and 3, and is aimed at giving a quantitative basis - as complete as possible - to model real 
search problems involving intermittent searchers. 

We first define our model and give general notations that we will use in this article. Then 
we systematically examine each case, studying the search problem in dimension 1, 2 and 3, 
where for each dimension different types of motion in the slow phase are considered. Each 
case is ended by a short summary, and we highlight the main results for each dimension. 
Eventually we synthesize the results in the table [Til where all cases, their differences and 
similarities are gathered. This table finally leads us to draw general conclusions. 
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II. MODEL AND NOTATIONS 



A. Model 

The general framework of the model relies on intermittent trajectories, which have been 
put forward for example in jo]. We consider a searcher that switches between two phases. 
The switching rate Ai (resp. A2) from phase 1 to phase 2 (resp. from phase 2 to phase 
1) is time-independent, which assumes that the searcher has no memory and implies an 
exponential distribution of durations of each phase i of mean = l/Aj. 

Phase 1 denotes the phase of slow motion, during which the target can be detected if it 
lies within a distance from the searcher which is smaller than a given detection radius a. a is 
the maximum distance within which the searcher can get information about target location. 
We propose 3 different modelings of this phase, in order to cover various real life situations. 

• In the first modeling of phase 1, hereafter referred to as the "static mode", the searcher 
is immobile, and detects the target with probability per unit time k if it lies at a 
distance less than a. 

• In the second modeling, called the "diffusive mode", the searcher performs a continuous 
diffusive motion, with diffusion coefficient D, and finds immediately the target if it 
lies at a distance less than a. 

• In the last modeling, called the "ballistic mode", the searcher moves ballistically in a 
random direction with constant speed Vi and reacts immediately with the target if it 
lies at a distance less than a. 

Some comments on these different model ings o f the slow phase 1 are to be made. The first 



two modes have already been introduced P, I2JI, while the analysis of the ballistic mode has 



never been performed. These 3 modes schematically cover experimental observations of the 



behavior of animals searching for food 
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2l| . where the slow phases of detection are often 



described as either static, random or with slow velocity. Several real situations are likely to 
involve a combination of two modes. For instance the motion of a reactive particle in a cell 
not bound to motors can be described by a combination of the diffusive and static modes. 
For the sake of simplicity, here we treat these modes independently, and our approach 
can therefore be considered as a limit of more realistic models. We note that a similar 
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version of the ballistic mode of detection has been discussed by Viswanathan et al. 
They introduced a model searcher performing randomly oriented ballistic movements (and 
of power law distributed duration), with detection capability all along the trajectory. For 
this 1 state searcher, the search time for a target (which is assumed to disappear after the 
first encounter) is minimized when the searcher performs a purely ballistic motion and never 
reorientates. In fact, the ballistic mode version of our model extends this model of 1 state 
ballistic searcher, by allowing the searcher to switch to a mode of faster motion, but with no 
perception. As it will be discussed, adding this possibility of intermittence enables a further 
minimization of the search time. Finally, combining these three schematic modes covers a 
wide range of possible motions, from subdiffusive (even static), diffusive, to superdiffusive 
(even ballistic). 

The phase 2 denotes the fast phase during which the target cannot be found. In this 
phase, the searcher performs a ballistic motion at constant speed V and random direction, 
redrawn each time the searcher enters a phase 2, independently of previous phases. In real 
examples correlations between successive phases could exist. If correlations are very high, 
it is close to a 1-dimensional problem with all the phases 2 in the same direction, a different 
problem already treated in |6||. We consider here the limit of low correlation, that is of a 
searcher with no memory skills. 

We assume that the searcher evolves in d-dimensional spherical domain of radius 6, with 
reflective boundaries and with one centered immobile target (bounds can be obtained in 



the case of mobile targets 



32, 



331). As the searcher does not initially know the target 



location, we start the walk from a random point of the d-dimensional sphere, and average 
the mean target detection time over the initial position. This geometry models the case 
of a single target in a finite domain, and also provides a good approximation of an infinite 
space with regularly spaced targets. Such regular array of targets corresponds to a mean- 
field approximation of random distributions of targets, which can be more realistic in some 
experimental situations. We note that in the 1 dimensional case, we have shown that a 
Poisson distribution of targets can lead to significantly different results from the regular 



distribution 



34 



35|. We expect this difference to be less in dimension 2 and 3, and we limit 



ourselves in this paper to the mean field treatment for the sake of simplicity. 
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B. Notations 



ti{r) denotes the mean first passage time on the target, for a searcher starting in the 
phase i from point r, where phase z = 1 is the slow motion phase with detection and phase 
2 = 2 is the fast motion phase without target detection. Note that in dimension 1, the 
space coordinate will be denoted by x, and in the case of a ballistic mode for phase 1, the 
upper index in tf stands for ballistic motion with direction ±x. The general method which 
will be used at length in this paper, consists in deriving and solving backward equations for 
ti{r) IsgI- These linear equations involve derivatives with respect to the starting position r. 

Assuming that the searcher starts in phase 1, the mean detection time for a target is then 
defined as : 

with Qd the d-dimensional sphere of radius b and V{fld) its volume. Unless specified, we 
will consider the low target density limit a <^h. 

Our general aim is to minimize as a function of the mean durations ri,r2 of each 
phase, and in particular to determine under which conditions an intermittent strategy (with 
finite T-i) is faster than a usual 1 state search in the phase 1 only, which is given by the 
limit Ti — > oo. In the static mode, intermittence is necessary for the searcher to move, 
and is therefore always favorable. In the diffusive mode, we will compare the mean search 
time with intermittence tm to the mean search time for a 1 state diffusive searcher t^i//, and 
define the gain as gain = tdiff/tm- Similarly in the ballistic mode, we will compare tm to the 
mean search time for a 1 state ballistic searcher thai and define the gain as gain = tf,ai/tm- 

Throughout the paper, the upper index "opt" is used to denote the value of a parameter 
or variable at the minimum of tm . 



III. DIMENSION 1 



Besides the fact that it involves more tractable calculations, the 1-dimensional case can 
also be interesting to model real search problems. At the microscopic scale, tubular struc- 



tures of cells such as axons or dendrites in neurons can be considered as 1-dimensional [25|. 
The active transport of reactive particles, which alternate diffusion phases and ballistic 
phases when bound to molecular motors, can be schematically captured by our model with 
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diffusive mode [2J|. At tfie macroscopic scale, one could cite animals like ants [37| which 
tend to follow tracks or one-dimensional boundaries. 



A. 



Static mode 



phase 2 




phase 1 



f 



V 
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X 



2a 



Figure 1: (Color online) Static mode in dimension 1 



In this section we assume that the detection phase is modeled by the static mode. Hence 
the searcher does not move during the reactive phase 1, and has a fixed reaction rate k per 
unit time with the target if it lies within its detection radius a. It is the limit of a very slow 
searcher in the reactive phase. 

1. Equations 

Outside the target (for x > a, we have the following backward equations for the mean 
first-passage time : 



Inside the target (x < a), the first two equations are identical, but the third one is written : 




(2) 




(3) 




(4) 




(5) 
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We introduce ^2 
equations : 



and to 



. Then outside the target we have the following 



dx T2 



-{t2-tl] 
Tl 



-1. 



Inside the target the first two equations are identical, but the last one writes : 



Tl 



k]ti = -1. 



(6) 
(7) 
(8) 

(9) 



Due to the symmetry x ^ — x, we can restrict the study to the part x G [0, a] and the part 
X G [a,b]. This symmetry also implies : 



dx 



dx 



x=0 



x=b 







0. 



In addition, continuity at a; = a for and ^2 gives: 

= a)= tl'^'ix = a) 



,d,in / \ id.out / \ 

12 [X = a) = 12 [X = a). 



(10) 
(11) 

(12) 
(13) 



This set of linear equations enables an explicit determination of ti, ^2, ^2 inside and outside 
the target. 



2. Results 

An exact analytical expression of the mean first passage time at the target is then given 
by : 

where 13 = a/ (kri)^^ + 1. 

In order to determine the optimal strategy, we need to simplify this expression, by ex- 
panding ( |T4l) in the regime b ^ a : 
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We make the further assumption <C 1 and obtain, using /3 > 1 : 



1 a2b ^2 

KTi 3V^T2 a 



(16) 



Since /5 > 1 and /? > l/(A;ri), we obtain in the Hmit h':^ a : 



(n + Ta) 



1 



1 



(17) 



This simple expression gives a very good and convinient approximation of the mean first 
passage time at the target as shown in FiglH 




r2/Tr ' 
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Figure 2: (Color online) Static mode in 1 dimension. Exact expression of (fT4l) (lines) compared 
to the approximation of (flTl) (symbols), both rescaled by tm* ([20]) . t°^* from (fTSj) . r2^* from 
(fT9l) . y = 1, A; = 1. h/a = 10 (green, squares), b/a = 100 (red, circles), h/a = 1000 (blue, crosses). 



We use this approximation (fTTl) to find ri and r2 values which minimize tr, 

1/4 



opt 



Vk \ 12a 



(18) 
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a b 
y V 3a' 



(19) 



It can be noticed than r2^* does not depend on k. Then the expression of the minimal value 



of the search time (flTl) with ti = t"^* and T2 



^opt 



IS : 



fopt 




2bk /3a^ 

w It 




(20) 



3. Summary 

For the static modeling of the detection phase in dimension 1, in the b ^ a limit, the 
mean detection time is : 



[Tl + T2) 



1 



(21) 



Intermittence is always favorable, and the optimal strategy is realized when r°^* = 
(lla) ^''^ and Tg^* = ^\J^- Importantly, the optimal duration of the relocation phase 
does not depend on A;, i.e. on the description of the detection phase. 



B. Diffusive mode 




phase 2 



2a 



Figure 3: (Color online) Diffusive mode in one dimension 



We now turn to the diffusive modeling of the detection phase. The detection phase 1 
is now diffusive, with immediate detection of the target if it is within a radius a from the 
searcher. 
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1. Equations 



Along the same lines, the backward equations for the mean first-passage time read outside 
the target {x > a) : 

^ d'^tl 1 ^ tr) 



and inside the target (x < a) : 

- ^4 = -1 (25) 

ax T2 

- V-^ 12 = -1 26) 

ti = 0. (27) 
We introduce the variables t2 = ^ind ^ = ^^-J^. This leads to the following system 



outside the target {x > a) : 



= ^4 (28) 

dx T2 



r2 dt2 , 1 



and inside the target (x < a) 



V'T2-^ + -ih-t2) = -l (29) 

dx T2 

„ d'^ti 1 , , , , 

= -ti„ (31) 

dx r2 

V^V^^ - -t2,.^ = -1 (32) 

fix T2 

ti = 0. (33) 

Interestingly, this system is exactly of the same type that what would be obtained with 2 
diffusive phases, with 02"^-^ = V'^T2 in phase 2. Boundary conditions result from continuity 
and symmetry : 

ti(a) = (34) 

tt{a) = tlM (35) 

^2 («) = (36) 
11 



dt2 
dx 

dti 
dx 

dt2,in 

dx 



x=b 



x=b 



x=0 



(37) 
(38) 
(39) 



2. Results 



Standard but lengthy calculations lead to an exact expression of mean first detection time 
of the target tm given in appendix IVII A 1[ We first studied numerically the minimum of tm 
in appendix IVII A 21 and identified 3 regimes. In the first regime {b < y ) intermittence is 
not favorable. For b > y intermittence is favorable and two regimes {-^^2 < 1 and > 1) 
should be distinguished . We now study analytically each of these regimes. 



3. Regime where intermittence is not favorable : b < y 

lib < y, the time spent to explore the search space is smaller in the diffusive phase than 
in the ballistic phase. Intermittence cannot be favorable in this regime, as confirmed by the 
numerical study in appendix IVII A2[ 

Without intermittence, the searcher only performs diffusive motion and the problem can 
be solved straightforwardly. The backward equations read tdj// = inside the target (x < a), 
and outside the target {x > a) : 

dx^ 

Since tdiff{x = a) = and = 0, we get U//(x) = j^{{b - a)^ - (6 - xf ). The 

x=b 

mean first passage time at the target then reads : 

f - - (41) 

U^n - (41) 

which in the limit b ^ a leads to : 

h2 



td^ff ^ (42) 
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4- Optimization in the first regime where intermittence is favorable : b < y and -^yi ^ 1 



6^2 



As explained in details in appendix IVII A 31 we use the approximation of low target 

density (6 ^ a), and we use assumptions on the dependence of r°^* and r2^* on b and a. 

These assumptions lead to the following approximation of the mean first passage time : 

h 1 



3VV| 



+ 



(43) 



We checked numerically that this expression gives a good approximation of in this regime, 
in particular around the optimum (FigHj). 




10"' 10" 10' 




Figure 4: (Color online) Diffusive mode in 1 dimension, from (j42]l . and tm exact 

expression (|158p (line), approximation in the regime of favorable intermittence and -^yi ^ 1 (l43l) 
(symbols), a = 1 and h = 100 (green, circles), a = 1, 6 = 10^ (red, crosses), a = 10, & = 10^ (blue, 
squares). D = 1, V = 1. r°^* is from expression ([441) . r2^* is from expression ((45]l 



The simplified expression (l43l) is minimized for : 

_ 1 ,[Wd 



^2 



9^4 



35 64 



This compares to the case without intermittence (|4T1) according to : 



gain 



opt 



tdiff 
.opt 



'24 fbV 



0.13 



6y 



These results are in agreement with numerical minimization of the exact tm (Table 
appendix IVII A 2^ . 



(44) 
(45) 
(46) 

(47) 
in 
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5. Optimization in the second regime where intermittence is favorable : b < y and 1 ^ 





Figure 5: (Color online) Diffusive mode in 1 dimension, -f^, tdiff from (j42]l . and tm exact 
expression (|158p (line), approximation in the regime of favorable intermittence and -^iyi *C 1 (j48l) 
(symbols), a = 10 and 6 = 100 (green, circles), a = 10, & = 1000 (red, crosses), a = 100, 6 = 10^ 
(blue, squares). D = 1, V = 1. r°^* is from expression (j49l) . r2^* is from expression ([50]) 



We start from the exact expression of (11581) . As detailed in appendix lVII A41 we make 
assumptions on the dependence of r°^* and r2^* with h and a, and use the assumptions that 
6 > a and 1 > It leads to : 



a ah 
[n + raj I , + ^775-3 



(48) 



This expression gives a good approximation of tm-, at least around the optimum (Figl5j), 
which is characterized by: 

Dh 

(49) 



^opt 



ASV^a 
a b 



- (^) 



3/2 



(50) 
(51) 

These results are in very good agreement with numerical data (TableHTTlin appendix lVII A 2^ . 
Note that the gain can be very large at low target density. 



gain 



1 aV b 
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6. Summ,ary 



We calculated explicitly the mean first passage time t,„ in the case where the detection 
phase is modeled by the diffusive mode. We minimized as a function of ti and r2, the 
mean phases durations, with the assumption a <^b. There are three regimes: 

• when 6 < y, intermittence is not favorable. r°^* oo, r2^* 0, t^* = tdijf — |^ 



when b> y and ^ 1, intermittence is favorable, with Tg^* = 2r°^* = -y/^?, and 

^opt 



• when b > y and <C 1, intermittence is favorable, with r°^* = 4^^^, ^2^* = -p^ y 3^, 

This last regime is of particular interest, since the value obtained for Tg^* is the same as in 
the static mode (cf. section [Til A 31) . 



C. Ballistic mode 




Figure 6: (Color online) Ballistic mode in one dimension 



We now treat the case where the detection phase 1 is modeled by the ballistic mode. 
This model schematically accounts for the general observation that speed often degrades 
perception abilities. Our model corresponds to the extreme case where only two modes are 
available : either the motion is slow and the target can be found, or the motion is fast and 
the target cannot be found. Note that this model can be compared with 
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1. Equations 



The backward equations read outside the target {x > a) : 



Vl 



dx 7-1^2 
(it If 



1 /tj , 



dx Ti \ 2 



1 ftt , tr 



-1. 



(53) 

(54) 
(55) 
(56) 



Defining t 



and ti 



v+t. 

2 """" ^ 2 

with Vl ,ti ^ t2, ^2 ti) ■ 



, we get the following equations (and similar expressions 



dti I , 
vi-^ + -{t2-ti) = -l 

dx Ti 

dti 1 w „ 
dx Ti 



which eventually lead to the following system : 

2 dHi , 1 



V,Ti 



+ -{t2-h) 



dx'^ Ti 



dX^ T2 



together with : 



dx 

Inside the target (x < a), one has t^'*"(x) = tf'*"(x) = 0, and : 



V 



dtn 



^2 

dx T2 
1 



-1 



dx T2 



-1. 



Finally, the boundary conditions read : 



dt2 
dx 



(57) 
(58) 

(59) 
(60) 

(61) 
(62) 

(63) 
(64) 

(65) 



x=b 
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dti 
dx 



x=b 



^2 («) = ttiJ 



a) 



t 



dt 



2, in 



dx 



= 



x=0 

t^{a) = 



(66) 
(67) 
(68) 
(69) 
(70) 



2. Results 



The exact expression of tm (cf. appendix IVII Bl) is obtained through lengthy but standard 
calculations. To simplify this expression, we consider the small density limit a/b and 
finally obtain the following very good approximation of tm (FiglT]) : 

(n + T2)b 



where 



a 



3/2 



(y/g - Li)(Li + L2) + y/a{L2 - Li + y^)X + ^^^2(^2 - 
((Li + v/^)X2 + (Li - v/^))(v/^+ L2 - Li) 

2a 

X = e^2 



/^2 = VT2 
Li = ViTi 





Vl = l 


•i;^ = 0.1 


vi = 0.01 


vi = 0.001 




6 = 5 


Ti ^ 00 


ri ^ 0, r"^* = 0.86 


^Opt,th ^ 


6 = 50 


Tl — > 00 


n ^ 0, r°^* = 2.9 


^opt,th ^ 2.9 


6 = 500 


Ti — > CX) 


n ^ 0, r°^* = 9.1 




6 = 5000 


Ti — > 00 


n ^ 0, rf * = 29 





(71) 

(72) 

(73) 
(74) 
(75) 
(76) 



Table I: Ballistic mode in 1 dimension. Numerical minimization of the exact tm (I200p . Values of ri 
and T2 at the minimum. Comparison with theoretical T2- a = 0.5, V = 1. 

A numerical analysis indicates (FigjTl Table [11) that, depending on the parameters, there 
are two possible optimal strategies : 
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1. X 10^ 




_opt 



T-2/T-2 

Figure 7: (Color online) Ballistic mode in 1 dimension. Comparison between low density approx- 
imation (fTTIl (symbols) and the exact expression of (I200p (line), (a), (b) : as a function of 



Ti, with T2 = O.lTg^* (red, crosses), T2 = t!^^'' (blue, squares), T2 = IOt!^^'' (green, circles), (c), (d) : 
tm as a function of T2/r2^*, with ri = (red, crosses), n = 10 (violet, circles), n = 100 (blue, 
squares), ri = 10000 (green, diamonds), (a), (c) : = 0.12 > vf : intermittence is not favorable, 
(b), (d) : vi = 0.06 < vf : intermittence is favorable. r2^* is from the analytical prediction ([8T]) . 
a = 1, y = 1, 6 = 100. 

• Ti ^ 00 . Intermittence is not favorable. 



opt 



opt 



Ti — > 0, r2 = Tg^*. Intermittence is favorable. 



We now study analytically these regimes. 
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3. Regime without intermittence ; ri — > oo 

In this regime, there is no intermittence. The searcher starts either inside the target {x 
in [—a, a]) and immediately finds the target, or it starts at a position x outside the target. 
We can therefore take x G [a,b]. If the searcher goes in the —x direction, it find its target 
after T = {x — a)/vi. If the searcher goes in the +x direction, it finds its target after 
T = {{b-x) + {b- a))/vi. This leads to : 

thai = T / dx = \ \ 77 

J a Vi bVi 

4- Intermittent regime 

We take the limit ri in the expression of tm ( TZTI) and obtain : 

}"V-=v[w.^:nz-J '''' 

Taking the derivative with respect to L2 yields : 

^ lim t„J oc 12ae^ + 2be^ - 6 - be^ (79) 



dL2 \n^0 
which has only one positive root : 



Lf = — (80) 

In M + 6a/b + 2^3a/b + da"^ /b^ 



In the limit fe ^ a it leads to : 



which is in agreement with numerical minimization of the exact mean detection time shown 
in the table [B 

The mean first passage time at the target is minimized in the intermittent regime for 
Ti ^ and T2 = r2^*. We replace T2 by (IHB in the expression (l78l) . and take 6 ^ a to finally 
obtain : 

gain=— (83) 

2 vi y b 

This shows that the gain is larger than 1 for vi < vf = ^-^\/f' '^^ich defines the regime 
where intermittence is favorable. 
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5. Summ,ary 



In the case where the phase 1 is modeled by the ballistic mode in dimension 1, we 
calculated the exact mean first passage time tm at the target, tm can be minimized as a 
function of Ti and r2, yielding two possible optimal strategies : 

• for vi > Vi = V^y/^, intermittence is not favorable : r"^* oo, r2^* 

• for Vi < vf = y^\f\-, intermittence is favorable, with r°^* and = -^sj^ 

Note that the model studied in [2|] shows that when targets are not revisitable, the 
optimal strategy for a 1 state searcher is to perform a straight ballistic motion. This strategy 
corresponds to ti ^ cxd in our model. Our results show that if a faster phase without 
detection is allowed, this straight line strategy can be outperformed. 

D. Conclusion in dimension 1 

Intermittent search strategies in dimension 1 share similar features for the static, diffusive 
and ballistic detection modes. In particular, all modes show regimes where intermittence is 
favorable and lead to a minimization of the search time. Strikingly, the optimal duration 
of the non-reactive relocation phase 2 is quite independent of the modeling of the reactive 
phase : r2^* = -^/sj^ for the static mode, for the ballistic mode (in the regime vi < vf ~ 

2"\/^)' diffusive mode (in the regime b > y and a ^ y\f^)- This shows the 

robustness of the optimal value r2^*. 



IV. DIMENSION 2 



The 2-dimensional problem is particularly well suited to model animal behaviors. It is 
also relevant to the microscopic scale, since it mimics for example the case of cellular traffic 
on membranes 25 1. The results for the static and diffusive modes, already treated in 0, ^, 
are summarized here for completeness. While in dimension 1 the mean search time can be 
calculated analytically, we introduce in dimension 2 (and later dimension 3) approximation 
schemes, which we check by numerical simulations. In these numerical simulations, diffusion 



was simulated using variable step lengths, as in 



38|, and we used square domains instead of 
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disks for numerical convenience (it was checked numerically that results are not affected as 
soon as 6 ^ a). 



A. Static mode 




Figure 8: (Color online) Static mode in dimension 2 



We study here the case where the detection phase is modeled by the static mode : the 
searcher does not move during the detection phase and has a finite reaction rate with the 
target if it is within its detection radius a. 



1. Equations 



The mean first passage time (MFPT) at a target satisfies the following backward equa- 



tions 



39|: 



27rri 



/ ihiV) - h{r))de^ - kiai^Mr) = -1. (84) 

Jo 

V ■ VM^) - -MT") - hiT")) = -1 (85) 

The function la is defined by la{~^) = 1 inside the target (if ll^] < a) and la(^) = 
outside the target (if |^| > a). In the present form, these integro-differential equations 
(completed with boundary conditions) do not seem to allow for an exact resolution with 
standard methods. ^2 is the mean first passage time on the target, starting from in phase 
2, with speed V , of angle 6^, and with projections on the axes Vx, Vy. i and j can take 
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either x ov y as a value. We use the following decoupling assumption 



{ViV,h)e, c:. {V^^,)e,{h)e^ (86) 

and finally obtain the following approximation of the mean search time, which can be checked 
by numerical simulations : 

tm = ^ {^(1 + kn)iy' - + \ [8y' + (1 + kn) (Vln(y/x) + (y' - - + 8))] } 

, (87) 

/ 2kTx a , / 2kTi h 
where x = 'i/ ; — -- — and?/='i/ ; — — — 88 

In that case, intermittence is trivially necessary to find the target: indeed, if the searcher 
does not move, the MFPT is infinite. In the regime 6 ^ a, the optimization of the search 
time (l87l) leads to : 

Tr = ^{Hhla)-ll2f'\ (90) 

and the minimum search time is given in the large volume limit by : 
52 2^/4 (a^ - 462) ln(6/a) + 2h^ - 



fopt _ 

a^k ^/Vk^ (21n(Va)- 1)3/4 



V2 (96a2 - 1926^) \n^{h/a) + (1926^ - lUa%'^) \n{h/a) + 460^62 _ 47^4 + ^4 

(91) 



A9,ahW (21n(6/a) - 1)3/2 

2. Summary 

In the case of a static detection mode in dimension 2, we obtained a simple approximate 
expression of the mean first passage time at the target. With the static detection mode, in- 
termittence is always favorable and leads to a single optimal intermittent strategy. The min- 
imal search time is realized for t"^^ = {y^Y^'^ ^ 2in(b/a)-i j ^opt _ a_ {\Yi{h/a) — 1/2)^^^. 
Importantly, the optimal duration of the relocation phase does not depend on k, i.e. on the 
description of the detection phase, as in dimension 1. 
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phase 2 




Figure 9: (Color online) Diffusive mode in dimension 2 



B. Diffusive mode 



We now assume that the searcher diffuses during the detection phase. For this process 



the mean first passage time at the target satisfies the following backward equation 

-2tv 



DVitAr 



1 



2'KT- 



{h{r)-t,{r))de^ = -i 



1 Jo 



V ■ vM^) - -ihiT") - = -1 

^2 



(92) 



(93) 



with ti(T^) = inside the target (r < a). We use the same decoupling assumption than for 
the static case (l86l) . It eventually leads to the following approximation of the mean search 
time : 



, 1 - ayt' 

[a Dti) 



with L± = Iq 



M L_ a'Dn (3-4 ln(6/a))6^ - 4a%^ + 



aaihya' - 1) 



2Li L. 



8Dto 



1,2 2 

— a 



(Ii(6a)Ki(aa) — Ii(aa)Ki(6a))±a y Dt2 I 



and M = lo 



Dt2, 



[\i{ha)Ko{aa) + Io(aa)Ki(6a))-4 



a'\jDr2 
a{b^ - a2)2 



Ii 



Dt2 



(94) 

[Ii(6a)Ko(aa;) + Io(aa)Ki(( 
[li{ba)Ki{aa) — Ii(aQ;)I 



where a = {1/{Dti) + 1/ (Dra))^/^ and D = v'^T2. We then minimize this time as a function 
of Ti and T2. 
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1. a <h <^ D/V : intermittence is not favorable 



In that regime, intermittence is not favorable. Indeed, the typical time required to explore 
the whole domain of radius b is of order b'^/D for a diffusive motion, which is shorter than 
the corresponding time b/V for a ballistic motion. As a consequence, it is never useful to 
interrupt the diffusive phases by mere relocating ballistic phases. We use standard methods 
to calculate the mean first passage time at the target in this optimal regime of diffusion 
only : 

D d f dt\ 

^^.]=-^ (95) 



r dr \ dr ^ 

The boundary conditions t{a) = and 4^(r = 6) = lead to 



t,,,f = [Aa%' -a'- 3b' + 46Mn ^ ) , (96) 



W/ = ^7^(-3 + 41n-) (97) 



and we find in the limit b ^ a : 

2. a <^ D/V <^ b : first regime of intermittence 

In this second regime, one can use the following approximate formula for the search time: 

= Jl^J'^J^ {41n(6/a) - 3 - 2^^{\^{aa) + 7 - ln2)} , (98) 

7 being the Euler constant. An approximate criterion to determine if intermittence is useful 
can be obtained by expanding in powers of 1 /ri when ri ^ oo (ri — oo corresponds to 
the absence of intermittence), and requiring that the coefficient of the term 1/ri is negative 
for all values of r^- Using this criterion, we find that intermittence is useful if 

72 exp(-7/4 + i)VblD - 4 ln(6/a) + 3 > 0, (99) 

In this regime, using Eq(l98l). the optimization of the search time leads to 



\? 4 In — 5 + c b -\/41nu; — 5 + c 



opt ^ " — - ' - opt ^ V ^ ^ ^ I . . 

^ Du;2(4in«;-7 + c)' ' V w ^ ' 



where w is the solution of the implicit equation w = 2Vbf{w)/D with 

-8(ln wf + (6 + 8 ln(6/a)) In w - 10 ln(6/a) + 11 - c(c/2 + 2 ln(a/&) - 3/2) 



V41nw — 5 + c 2 



(101) 
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and c = 4(7 — ln(2)), 7 being the Euler constant. A useful approximation for w is given by 

2Vb / Vb \ , , 



D ' \2D\n{b/a 
The gain for this optimal strategy reads : 

tdiff 1 41n6/a-3 + 4aV62 _a4/^4 / ^ ^jj 
gain = -if ^ TTTT-TTT ^^^TTTT— TT-TTT— T TT— r + ' — ' 



2 41n6/a-3 + 2(41nu;)ln(6/au7) V41ntf;-5 bV {Alnw - 5^/^ J 

(103) 

If intermittence significantly speeds up the search in this regime (typically by a factor 2), it 
does not change the order of magnitude of the search time. 



3. D/V <^ a <^b : "universal" regime of intermittence 

In the last regime D/V ^ a <^ b, the optimal strategy is obtained for 

and the gain reads: 

t.ff^_ V2av ( 1 Io(2/v/21n(Va)-i; ^ 

8D 1^41n(6/a)-3i^ (2/ ^2 ln(6/a) - l) 2^2\n{b/a) - 1 ^ 

(105) 

Here, the optimal strategy leads to a significant decrease of the search time which can be 
rendered arbitrarily smaller than the search time in absence of intermittence. 



-1 



4- Summary 

We studied the case where the detection phase 1 is modeled by a diffusive mode, and 
obtained an approximation of the mean first passage time at the target. We found that 
intermittence is favorable (i.e. better than diffusion alone), in the regime of large system 
size b ^ D/V. The optimal intermittent strategy then follows two subregimes : 

• if a ^ D/V, the best strategy is given by (11 001) . The search is significantly reduced 
by intermittence but keeps the same order of magnitude as in the case of a 1-state 
diffusive search. 

• if a D/V, the best strategy is given by (11040 . and weakly depends on b. In this 
regime, intermittence is very efficient as shown by the large gain obtained for V large. 
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c. 



Ballistic mode 



^ phase 1 



'\ phase 2 



y 




X 



V, 



Figure 10: (Color online) Ballistic mode in two dimensions 



In this case, the searcher has access to two different speeds: one (V) is fast but prevents 
target detection, and the other one (vi) is slower but enables target detection. 

1. Simulations 

Since an explicit expression of the mean search time is not available, a numerical study 
is performed. Exploring the parameter space numerically enables to identify the regimes 
where the mean search time is minimized. Then, for each regime, approximation schemes 
are developed to provide analytical expression of the mean search time. 

The numerical results presented in figure [TT] suggest two regimes defined according to a 
threshold value of vi to be determined later on : 

• for Vi > vf, tm is minimized for T2 ^ 

• for Vi < ff, tm is minimized for Ti — > . 

2. Regime without intermittence (t2 — > 0, ri — > ooj 

Qualitatively, it is rather intuitive that for vi large enough (the precise threshold value vf 
will be determined next), phase 2 is inefficient since it does not allow for target detection. 
The optimal strategy is therefore r2 in this case. In this regime, the searcher performs 
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Figure lli (Color online) Ballistic mode in 2 dimensions. In(tj7^ ) as a function of ln(r2). Simulations 
(symbols), diffusive/diffusive approximation jM]) with (|109p (colored lines), ri — > limit (|110p 
(black line), ti ^ oo (no intermittence) (|108p (dotted black line). 6 = 30, a = l,y = l. ti = 
(black, stars), ri = 0.17 (yellow, solid circles), ti = 0.92 (green, diamonds), ri = 5.0 (blue, X), 
Ti = 28 (purple, circles), ri = 150 (red, +), ri = 820 (brown, squares). 
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Figure 12: (Color online) Ballistic mode in 2 dimension, tm* as a function of b, logarithmic scale. 
Regime without intermittence {t2 = and n — > oo, = 1), analytical approximation (|108|) 
(blue line), numerical simulations (blue circles). Regime with intermittence (with ri = 0, F = 1), 
analytical approximation (11121) (red line), numerical simulations (red squares), a = 1 



a ballistic motion, which is randomly reoriented with frequency I/ti. Along the same line 
as in [2| (where however the times between successive reorientations are Levy distributed), 
it can be shown that the optimal strategy to find a target (which is assumed to disappear 
after the first encounter) is to minimize oversampling and therefore to perform a purely 
ballistic motion. In our case this means that in this regime T2 ^ 0, the optimal ti is given 
by T^^'^ 00. 
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In this regime, we can propose an estimate of the optimal search time tbai- The surface 
scanned during 6t is 2avi6t. p{t) is the proportion of the total area which has not yet been 
scanned at t. If we neglect correlations in the trajectory, one has : 

dp 2avip{t) 



dt nb^ 
Then, given that p{t = 0) = 1, we find 



(106) 



p{t) = exp ( ) , (107) 



TTb^ 

and the mean first passage time at the target in these conditions is : 

r dp , vr62 

tbal = - / t-^dt = . 108 

Jo dt 2avi 

This expression yields (Fig JT2l) a good agreement with numerical simulations. Note in par- 
ticular that tbal oc 



3. Regime with intermittence ti ^ 

In this regime where vi < vf, the numerical study shows that the search time is minimized 
for Ti ^ (FigJTTl). We here determine the optimal value of r2 in this regime. To proceed 
we approximate the problem by the case of a diffusive mode previously studied (l94l) . with 
an effective diffusion coefficient : 

D='-ti. (109) 
This approximation is very satisfactory in the regime ri — > as shown in FigJTIl 



We can then use the results of 



ictorv 

0,3 

_f a^\( l (3 + 41n(f))6^-4a^6^ + a^ a ( .VH^ 



in the ri ^ regime and obtain: 

4 . /U2 \ iJ^ 



The calculation of r2^* minimizing tm then gives: 




(110) 



(111) 
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In turn, replacing T2 by Tg^* (11 111) in (II 101) . we obtain the minimal mean time of target 
detection : 



,ort^^(,_-l\ L .. (3 + 41n(f))fe^-4a^&^ + a^ n (^^ - 1 j Jp (2n) 
uV2V V &V I (62 - a2) ^ Ji (2n) 

'(112) 

with u = (in (^) — 1) ^ . It can be noticed that t^* oc Note that if 6 ^ a this last 
expression can be greatly simplified: 



aV y \a 

Finally the gain reads (using (11081) and (11131) ): 



C'^^Jlnf^^) (113) 



Numerical simulations of FigJT2l shows the validity of these approximations. 



4- Determination of vf 

It is straightforward than < V. Indeed, ii vi = V , phase 2 is useless, since it prevents 
target detection. Actually, an estimate of vf can be obtained from (1114^ as the value of vi 
for which gain = 1 : 

ttV / V 



We note that this expression (FigJT3l) gives the correct dependence on b, but it however 
departs from the value obtained by numerical simulations. This is due to fact that the 
expression of t^* with intermittence (I113P is an under estimate, while t^ai given in (llOSp is 
an upper estimate. It is noteworthy that intermittence is less favorable with increasing b. 
This effect is similar to the 1 dimensional case, even though it is less important here. It can 
be understood as follows: at very large scales the intermittent trajectory is reoriented many 
times and therefore scales as diffusion, which is less favorable than the non intermittent 
ballistic motion. 
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Figure 13: (Color online) Ballistic mode in two dimension, as a function of ln(6) by simulations 
(symbols), predicted expression (jllSp (red dotted line), predicted expression multiplied with a fitted 
numerical constant (blue line). V = 1, a = 1. 

5. Summary 

We studied the case of ballistic mode in the detection phase 1 in dimension 2. When 
vi > vf, the optimal strategy is to remain in phase 1 and to explore the domain in a purely 
ballistic way. Therefore, Tg^* , t"^*^ oo. When vi < vf, we find on the contrary 
rf^* — > and To^* = tt \ / In (-) — The threshold value is given by vf oc , ^ and shows 
that when the target density decreases, intermittence is less favorable. 

D. Conclusion of the 2-dimensional problem 

Remarkably, for the three different modes of detection (static, diffusive and ballistic), 
we find a regime where intermittence permits to minimize the search time for one and the 



intermittent strategies are robust and widely independent of the details of the description 
of the detection mechanism. 

V. DIMENSION 3 

The 3 dimensional case is also relevant to biology. At the microscopic scale, it corresponds 
for example to intracellular traffic in the bulk cytoplasm of cells, or at larger scales to animals 



same r2 , given 




As in dimension 1, this indicates that optimal 
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4l|. As 



living in 3 dimensions, such as plankton j4o| . or C.elegans in its natural habitat 
it was the case in dimension 2, different assumptions have to be made to obtain analytical 
expressions of the search time. We checked the validity of our assumptions with numerical 
simulations using the same algorithms as in dimension 2. 



A. Static mode 




Figure 14: (Color online) Static mode in three dimensions 



We study in this section the case where the detection phase is modeled by the static 
mode, for which the searcher does not move during the detection phase and has a finite 
reaction rate with the target if it is within a detection radius a. 



1. Equations 

Denoting tiir) the mean first passage time at the target starting from a distance r from the 
target in phase 1 (detection phase), and t2fi,^{r) the mean first passage time at the target 
starting from a distance r from the target in phase 2 (relocation phase) with a direction 
characterized by 6 and 0, we get : 

V^t2,e,^ + -{ti-t2,e,^) = -l- (116) 

T2 

Then one has outside the target (r > a) : 

— {-^ I ddsind ! d(t)t2fi.^-t^ = -I, (117) 



Ti V47r 



Jo 
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and inside the target (r < a) : 



1 1 r z"^" /I 

— — / dOsinO / d(j)t2e<h — ( \- k ] ti = —1. 

Ti^TTJo Jo ' ' \n 



(118) 



With t2 = ^ dOsine /q rf0t2,e,<A, 

we obtain outside the target (r > 
1 



{t2 - h) 



(119) 



and inside the target (r < a) 



-t2 - (- + k]ti = -1. 



(120) 



Making a similar decouphng approximation as in dimension 2, we finally get : 

-1. 



^At2--{h-t2) 



(121) 



We solve these equations for inside and outside the target, using the following boundary 
conditions : 

Jj-OUt 

2 = (122) 



dr 



dr 



and the condition that f^{0) should be finite. 



2. Results 



dr 



(123) 
(124) 



We find an explicit expression of the mean search time 
1 



tm = in + T2) 



with 



kn b^v 



L^(-2 63(6^_a^) + (53_,3) ^3 



P 



— sinh (a) + a cosh (a) b'^ 
a (— sinh (a) + a cosh (a)) 



and a 



3t^ 



In the limit b ^ a, this can be simplified to : 



+ /3] +-(6'-a') 



tm. = (n + r2) I h 



— sinh (a) a'^ + a cosh (a) 6'^ 9^2 



kri r|V2 V (~ ^i'^^ ('^) + ^^^^^ i^)) ^ 



a" 



(125) 
(126) 
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Assuming further that a is small, we use the expansion /3 ~ ^ (1 -tanh(a)/a)"^ ~ 
— (-^ + I) and rewrite mean search time as : 

a V 5 / 



This expression of tm can be minimized for 



^opt 



3 \ * la 



10 J \ Vk 

a 



and the minimum mean search time reads finally: 



3. Comparisons with simulations 



(129) 



.opt ^ ^12^, (130) 



.opt _ f^2A^'^ + ri3n 




Figure 15: (Color online) Static mode in 3 dimensions. In(tm) as a function of ln(r2) for different 
values of Ti, a and h/a. Comparison between simulations (symbols), analytical expression (|125|) 
(line), expression for 6 » a (11271) (small dots), and simple expression for 6 o and a small (|128|) 
(dotted line), n ~ r^^* ~ 0-74-^^ (USHl) (blue, crosses), n = 0.25^^ (red, circles), n = 2.5^^ 
(green, squares). V = 1, k = 1. 
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Data obtained by numerical simulations (FigHH and additionally in appendix IVII CI) are 
in good agreement with the analytical expression (|125p . In particular, the position of the 
minimum is very well approximated, and the error on the value of the mean search time at 
the minimum is close to 10%. Note that the very simple expression (11281) fits also rather 
well the numerical data, except for small r2 or small b. 

4- Summary 

In the case of a static detection mode in dimension 3, we obtained a simple approximate 
expression of the mean first passage time at the target tm = ^ '''^1"^'^^'' + 5t^v^ ) • 

has a single minimum for r°^* = (^) y^y^ ^md r2^* = \/h2y, and the minimal mean search 
time is -^1^ ^y^24^/^ + S^/^'j . With the static detection mode, intermittence is always 
favorable and leads to a single optimal intermittent strategy. As in dimension 1 and 2, the 
optimal duration of the relocation phase does not depend on k, i.e. on the description of the 
detection phase. In addition, this optimal strategy does not depend on the typical distance 
between targets b. 

One can notice than for the static mode in the three cases studied (1,2, and 3 dimensions). 



we have the relation : = y / {2k). The optimal durations of the two phases are related 
independently of the dimension. 

B. Diffusive mode 





z 



Figure 16: (Color online) Diffusive mode in three dimensions 
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We now study the case where the detection phase is modeled by a diffusive mode. During 
the detection phase, the searcher diffuses and detects the target as soon as their respective 
distance is less than a. 



1. Equations 



One has outside the target (r > a) : 



1 



T2 



D /\ti + ^ dOsine c/0t2,e,<^ - ^i^ = -1 

and inside the target (r < a) : 

1 



-t 
ti = 0. 



2,0,0 



With t2 = n dOsinO J^"" d(f)t2,e,<t), we get outside the target (r > a) 



D A * + — (t^* - *) = -1 



(132) 
(133) 

(134) 
(135) 

(136) 



The decoupling approximation described in previous sections then yields outside the target 



1/2 1 

— - A + — (t°"* - t°"*) 



r2 



(137) 



and inside the target (r < a) : 



V^T2 



Af"* 



Ant -| 

2 '-2 ~ 

T2 



(138) 



These equations are completed by the following boundary conditions : 

= 



dr 



dt"^* 
dr 



dtf^ 
dr 



(139) 

(140) 
(141) 
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2. Results in the general case 



Through standard but lengthy calculations we can solve the above system and get an 
analytical approximation of tm (cf. appendix IVII D if l. In the regime b ^ a, we use the 
assumption ^ (tiD)^^ + 3{t2v)^'^ <^ b and obtain : 

b^4{Ti + T2) tanh(K2a) + ^ 



^1 KiK^TiDa ^tanh(/€2a,) + — tanh(fi;2a) 



(142) 



with Ki = cind K2 = y^- As shown in FigHTIleft or in the additional Figl24l in 

appendix lVII D 21 tm only weakly depends on ri, which indicates that this variable will be less 
important than T2 in the minimization of the search time. The relevant order of magnitude 
for r°^* can be evaluated by comparing the typical diffusion length i^di/f = y/ODt and the 
typical ballistic length L^ai = Vt. An estimate of the optimal time r°*^* can be given by the 
time scale for which those lengths are of same order, which gives : 

ft n 

~ ^. (143) 

Note that taking ri = does not change significantly t^* (FigHTI left), and permits to 
significantly simplify tm '■ 



— tanh 

VHi \Vt2 \Vt2 



^"^ = ^l^-tanh|^l 1 (144) 



In turn, the minimization of this expression leads to : 

opt 



-\/3a 



(145) 



with X solution of : 



This finally yields : 



2tanh(a;) - 2a; + xtanh(a;)^ = 0. (146) 



r"^* ^ 1.078:^. (147) 

Importantly this approximate expression is very close to the expression obtained for the 
static mode (rg^* = \J\^ — 1.095-^) (11301) . and there is no dependence with the typical 
distance between targets b. The simplified expression of the minimal tm (jl44l) can then be 
obtained as: 

C* = J^J— {x - tanh(x))"^ ^ S.IS^TT (148) 
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and the gain reads: 



9a^n = ^ ^ 0.15^. (149) 



3. Comparison between analytical approximations and numerical simulations 




b/a b/a 
Figure 17: (Color online) Diffusive mode in 3 dimension. Comparison between analytical approxi- 
mations p43|) (|145|1 (black lines) and numerical simulations : the symbols are the values of n and 

T2 for which < l_Qr^^opt,s^mulat^on ^ ^ ^ ^qq^ V = l,D=l 

Numerical simulations reveal that the minimum of with respect to Ti is shallow as it 
was expected (cf FigHT] left) . It approximately ranges from to the theoretical estimate 
f ll43p . The value 7-°^'*''**'" at the minimum is close to the expected values (1145^ (cf. FigJTTi 
right), except for very small 6, which is consistent with our assumption 6 ^ a. We can 
then conclude than the position of the optimum in ri and T2 is very well described by the 
analytical approximations, even if the value of at the minimum is underestimated by our 
analytical approximation by about 10-20% (FigJTSll. 



^. Case without intermittence : 1 state diffusive searcher 



If the searcher always remains in the diffusive mode, it is straightforward to obtain (cf. 
appendix [VTins]): 

tdiff = TFTTT^ + 56' - 96^a - a!") , (150) 



which gives in the limit 6/a ^ 1 
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Figure 18: (Color online) Diffusive mode in 3 dimension, tm/tdiff (tdiff given by (|150p ) as a 
function of T2 for different values of the ratio b/a (logarithmic scale). The full analytical form 
(|238p (plain lines) is plotted against the simplified expression (|142p (dotted lines), the simplified 
expression with ri = (11441) (small dots), and numerical simulations (symbols) for the following 
values of the parameters (arbitrary units): a = 1 (green, squares), a = 5 (blue, stars), a = 7 (purple, 
circles), a = 10 (red, +), a = 14 (brown, X), a = 20 (orange, diamonds), ri = 6 everywhere except 
for the small dots, v = 1, D = 1. tm/tdiff presents a minimum only for o > Cc ^ 4. 

5. Criterion for intermittence to he favorable 

There is a range of parameters for which intermittence is favorable, as indicated by FigJTSl 
Both the analytical expression for t"^* in the regime without intermittence (11501) and with 
intermittence (1148P scale as It'. However, the dependence on a is different (cf. appendix 
IVIID4D . In the diffusive regime, tm oc a , whereas in the intermittent regime tm oc a . 
This enables to define a critical Oc, such that when a > ac, intermittence is favorable: 
ttc — 6.5y is the value for which the gain (|149p is 1. 

6. Summary 

We studied the case where the detection phase 1 is modeled by a diffusive mode, and 
calculated explicitly an approximation of the mean first passage time at the target. We 
found that intermittence is favorable (i.e. better than diffusion alone) when a > ac — 6.5y : 
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• if a < Oc, the best strategy is a 1 state diffusion, without intermittence, and the mean 
first passage time at the target is tm — 375^- 

• if a > ac, intermittence is favorable. The dependence on ti is not crucial, as long as 
it is smaller than The value of T2 at the optimum is Tg^* ~ 1.08-^. The minimum 
search time is then t^* ~ 2.18^. 




Figure 19: (Color online) Ballistic mode in dimension 3 



We now discuss the last case, where the detection phase 1 is modeled by a ballistic 
mode. Since an explicit analytical determination of the search time seems out of reach, we 
proceed as in dimension 2 and first explore numerically the parameter space to identify the 
regimes where the search time can be minimized. We then develop approximation schemes 
in each regime to obtain analytical expressions (more details are given in appendix IVII Ep . 

1. Numerical study 

The numerical analysis puts forward two strategies to minimize the search time, depend- 
ing on a critical value to be determined (Fig l20l) : 

• when vi > vf, r°^* 00 and t^* — > 0. In this regime intermittence is not favorable. 

• when vi < vf, r°^* 0, and Tg^* finite. In this regime the optimal strategy is inter- 
mittent. 
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Figure 20: (Color online) Ballistic mode in 3 dimensions, tm as a function of T2 in loglogscale. 
Simulations (symbols). Approximation viTi < a (|238l with I253P (colored lines), approximation 
Ti = (I153|) (black line), approximation ti = and 6 S> a (11441) (dotted black line). Ballistic limit 
(r2 and ri ^ 00) (no intermittence) (I152p (green dotted line). (a),(d) : > vf, ri^i = 0.04, 
ri,2 = 0.2, ri,3 = 1, ri,4 = 5, ri,5 = 25. (b),(e) : ~ vf, ri,i = 0.08, ri,2 = 0.4, ri,3 = 2, ri,4 = 10, 
ri,5 = 50. (c),(f) : vi < vf, ri,i = 0.2, ri,2 = 1, n^3 = 5, ri,4 = 25, ri,5 = 125. F = 1, a = 1. n = 
(black, diamond), ri = rij (brown, +), ri = ri^2 (red, squares), ri = ri^3 (pink, stars), ri = ti^4 
(blue, circles), ti = ri^s (green, X) 

^. Regime without intermittence ( 1 state ballistic searcher) ; r2 — > 



Following the same argument as in dimension 2, without intermittence the best strategy 
is obtained in the limit ri — > cxd (cf . appendix IVII E Ij) in order to minimize oversampling of 
the search space. Following the derivation of 11081 (see appendix for details), it is found that 
the search time reads : 



46^ 

'ia?vi 



(152) 



3. Regime with intermittence 



In the regime when intermittence is favorable, the numerical study suggests that the best 
strategy is realized for ri ^ (Figl20l). In this regime ri ^ 0, the phase 1 can be well 
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approximated by a diffusion witli effective diffusion coefficient Deff (see (12530 ). We can 
tlien make use of tfie analytical expression tm derived in (|238l) . We therefore take ri = in 
the expression of tm (12381) , which yields : 

t^^r, = 0) ^ f ^ (5b'a' - 3b' - 2a') + ^^^'""'^'^ ] , (153) 

^ ' b^aV > V3a(n - tanh(n)) ] ^ ' 

where u = In the limit b ^ a, this expression can be further simplified (see (11441) ) to : 



— z — ^ tanh 



and one finds straightforwardly that Tg^* = where x is solution of a;tanh(x)^ + 

2tanh(x) — 2x = 0, that is x ~ 1.606. Using this optimal value of r2 in the expression 
of tm P44l) . we finally get : 

" 7!tanh(x)2^ - ^-^^^^ ^^^^^ 
These expressions show a good agreement with numerical simulations (Figl20 l Fig l2TI) . 



4- Discussion of the critical value vf 
The gain is given by : 

^a^^= -0-6177- (156) 

tm Vl 

As in dimension 2, it is trivial that vf < V, and the critical value vf can be defined as the 
value of vi such that gain = 1. This yields 

vf ~ 0.6V". (157) 

Importantly, vf neither depends on b nor a. Simulations are in good agreement with this 
result (cf. appendix IVIIE2I) . except for a small numerical shift. 



5. Summary 

We studied the case where the detection phase 1 is modeled by a ballistic mode in 
dimension 3. We have shown by numerical simulations that there are two possible optimal 
regimes, that we have then studied analytically : 
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Figure 21: (Color online) Ballistic mode in 3 dimension, tm as a function of 6, logarithmic scale. 
Regime without intermittence {t2 = and ri — > oo, vi = 1), analytical approximation (|152p 
(blue line), numerical simulations (blue circles). Regime with intermittence (with ri = 0, F = 1), 
analytical approximation (1155P (red line), numerical simulations (red squares), a = 1 



• in the first regime Vi > vf, the optimal strategy is a 1 state ballistic search (ri ^ oo, 
T2 = 0) and tm^Su; 

• in the second regime vi < vf, the optimal strategy is intermittent (ri = 0, T2 ~ l-l-^), 
and tm — 2.18^ (in the limit b 3> a). 

The critical speed is obtained numerically as vf ~ 0.5^ (analytical prediction : vf ~ 0.61^). 
It is noteworthy that when b ^ a, the values of Ti and T2 at the optimum, and the value of 
vf do not depend on the typical distance between targets b. 



D. Conclusion in dimension 3 

We found that for the three possible modelings of the detection mode (static, diffusive 
and ballistic) in dimension 3, there is a regime where the optimal strategy is intermittent. 
Remarkably, and as was the case in dimension 1 and 2, the optimal time to spend in the 
fast non-reactive phase 2 is independent of the modeling of the detection mode and reads 
Tg^* ^ l.ly. Additionally, while the mean first passage time on the target scales as b^, the 
optimal values of the durations of the two phases do not depend on the target density a/b. 
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VI. DISCUSSION AND CONCLUSION 



The starting point of this paper was the observation that intermittent trajectories are 
observed in various biological examples of search behaviors, going from the microscopic 
scale, where searchers can be molecules looking for reactants, to the macroscopic scale of 
foraging animals. We addressed the general question of determining whether such kind 
of intermittent trajectories could be favorable from a purely kinetic point of view, that is 
whether they could allow to minimize the search time for a target. On very general grounds, 
we proposed a minimal model of search strategy based on intermittence, where the searcher 
switches between two phases, one slow where detection is possible, the other one faster but 
preventing target detection. We studied this minimal model in dimensions 1, 2 and 3, and 
under several modeling hypotheses. We believe that this systematic analysis can be used as 
a basis to study quantitatively various real search problems involving intermittent behaviors. 

More precisely, we calculated the mean first passage time at the target for an intermittent 
searcher, and minimized this search time as a function of the mean duration of each of the two 
phases. The table HTl summarizes the results. In particular, this study shows that for certain 
ranges of the parameters which we determined, the optimal search strategy is intermittent. 
In other words, there is an optimal way for the intermittent searcher to tune the mean time 
it spends in each of the two phases. We found that the optimal durations of the two phases 
and the gain of intermittent search (as compared to 1 state search) do depend on the target 
density in dimension 1. In particular, the gain can be very high at low target concentration. 
Interestingly, this dependence is smaller in dimension 2, and vanishes in dimension 3. The 
fact that intermittent search is more advantageous in low dimensions (1 and 2) can be 
understood as follows. At large scale, the intermittent searcher of our model performs 
effectively a random walk, and therefore scans a space of dimension 2. In an environment 
of dimension 1 (and critically of dimension 2), the searcher therefore oversamples the space, 
and it is favorable to perform large jumps to go to previously unexplored areas. On the 
contrary, in dimension 3, the random walk is transient, and the searcher on average always 
scans previously unexplored areas, which makes large jumps less beneficial. 

Additionally, our results show that, for various modeling choices of the slow reactive phase, 
there is one and the same optimal duration of the fast non reactive phase, which depends 
only on the space dimension. This further supports the robustness of optimal intermittent 
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search strategies. Such robustness and efficiency - and optimality - could explain why 
intermittent trajectories are observed so often, and in various forms. 
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VII. APPENDIX 

A. Diffusive mode in dimension 1 

1. Exact results (cf. section UlIB ^|) 

1 in + T2) Num 

'- = -3 (3% Den ^'^^^ 

With : 

Num = «! + 02 + as + Q!4 + as + ae + 0:7 (159) 

Den = 71 + 72 + 73 + 74 (160) 

ai = (^(SLs' [Ll - Ll) + 2h^p) h^/pS + 3L1L2 (L2' - 2h^p) (161) 

a2 = -LihLl {2(3 + 3Ll) RC (162) 

as = Ll {2h^(3^ - 3Ll) BC (163) 

ai = h^^/p[QLl + h^p[l3 + Ll))RS (164) 

= ^hLl {Ah^(3 + 3Ll [lI - L^) BS (165) 

as = LiLl (3 {2h'^L2(3 + L^ B + h {3Ll {(3 + Ll) + 2h^/3) R) (166) 

ar = -Ll {3Ll + 2h^(3^) (167) 

'yi = LlLiR{C -1) (168) 

72 = v^/i {2Ll + Ll) RS (169) 

Js = ^/PLUB-1)S (170) 

74 = 2h(3Li {BC - 1) (171) 
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B = cosh ( ^"j (172) 



C = cosh [ 2h\/ + L^"^ ] (173) 



R = smh(^] (174) 



S = sinh \^2h^JL^^ + L^^ j (175) 
/5 = L? + L2 (176) 
Li = /D^ (177) 
L2 = Vt2 (178) 
h = b-a (179) 

^. Numerical study (cf. section UlIB ^) 

We studied numerically the optimum of the exact tm expression (|158p (Table UTTl) . We 
could distinguish 3 regimes : one with no intermittence, and two with favorable intermit- 
tence, but with different scalings. Intermittence is favorable when b > y. The demarcation 
line between the two intermittent regimes is -^^2 = 1- 

3. Details of the optimization of the regime where intermittence is favorable, with -^iyi ^ 1 



(cf. section \III B 4l ) 



We suppose that target density is low : | -C 1. 

We are interested in the regime where intermittence is favorable. We have both 2(6 — 
a)^/LY'^ + L^^ > 2^ and 2(6 - a)^/Lf^~^^L^ > 2^. In a regime of intermittence, one 
diffusion phase does not explore a significant part of the system : 6/Li ^ 1. Alternatively, 
having a ballistic phase of the size of the system is a waste of time, thus close to the optimum 
6/L2 > 1. Consequently 2(6 - a)^/L^'^ + L^'^ > 1. 

We use the numerical results (Table IIIII ) to make assumptions on the dependence of r°^* 
and Tg^* with the parameters. We define ki and k2 : 



ri = {ki)-'{yr) (180) 
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29000 


91000 


2.9.10^ 


9.2.10^ 


2.9.101^ 


9.8.10'^ 




£/i,2 
^2 


29000 


91000 


2.9.10^ 


9.1.10^^ 


2.9.10"^ 


9.1.10'= 



Table III: Diffusive mode in 1 dimension. Optimization of tm as a function of ri and T2 for different 
sets of parameters (D = 1, V = 1). For each (a, 6), numerical values for the exact analytical 
function pSSp are given with the values expected in the regimes where intermittence is favorable, 
either with > 1 {th,l), or with < 1 {th,2). gain^^'^ (j47|, gain = t^^/tdifj, gain^'*^ 
(1521). rf dm), r°^*, rf (09]). rf gSD, r"^*, rf Colors indicate the regime : red when 
intermittence is not favorable, green in the ^> 1 regime, blue in the -^/^ *C 1 regime. 

r2 = (fc2)-^(^)' (181) 
We make a development of for a. We suppose ki and /c2 do not depend on 6/a : 

We checked that this expression gives a good approximation of tm in this regime, in 
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particular around the optimum (FigHl in section [Til B 41) . 

Derivatives of ({431) as a function of ki and k2 must be equal to at the optimum. It leads 
to : 

-3kf +3k^ + 3^ik2 = (183) 
3kl - 2kl - kikl = (184) 
On four pairs of solutions, only one is strictly positive : 



opt 1 s^bW 



-r={/^ (186) 



4- Details of the optimization of the universal intermittent regime -^jTyz ^ 1 (cf- section UlI B 5\) 

We start from the exact expression of tm (11581) . We have to make assumption on the 
dependency of r2^* with b and a. We define / by r2 = j^sj^, and we suppose that / is 
independent from a/h. We make a development of a/6 — > 0. The first two terms give : 

b f 1 ^ \a + apWm[f 

^1 \ — 7^ (187) 



a\y 3Vf M a + 



This expression gives a very good approximation of in the hD"^ / {a^v"^) ^ 1 regime, 
especially close to the optimum (FiglHl in section [Til B 51) . 

We then minimize tm ( 1187^ as a function of / and ri. We introduce w defined as : 

aV fa 

» = -^- (188) 

We make an assumption on the dependency of r°^* with a/b, inferred via the numerical 
results : 

1 D b , , 

s = - — - 189 
Ti V'' a 

We write the equation p87p with these quantities. Its derivatives with / and s should equal 
zero at the optimum. It leads to : 

- v^s'^/V + V3s^/^-w^f^ + V3s wf^ + 6 y^wf^ + Qf^ = (190) 
6 s^'^w'^f + 6 s^/'^f n? - n?s^^ + 3 w;s / + 12 w;s + 6 ^sf = (191) 
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We take the equation (I190p and here we need to make the assumption than : a <C & <^ 



We get : 

V^s^/^w\f - 1) = 

Consequently / = 1. We incorporate this result to equation (11911) : 

12 - w^s^ + Ihws + 6 = 



(192) 



(193) 



The relevant solution is : 



^sol 



'l ^ 5v^+ 16m; 4 



3 w 



+ 



u 



With : 



u = i 270W + 27V3 + 192 w^Vs + 9 V 55 VSw + 84 u;2 + 27 



w 



(194) 



(195) 



When w ^ oo, Ssoi = 48. As we made the assumption = w ^ ^ 1, difference from the 
asymptote will be small (Figl22l). 



^sol 

64- 

62- 

60- 

58- 

56 

54 

52 

50 



Figure 22: (Color online) Diffusive mode in 1 dimension. Ssol (11941) (red line) as a function of ln(?i;), 
with the asymptote (blue dotted line) 



It leads to 



Tr 



opt 



lab 



a 




2 fopt V V 3a Vy 3a 
Db Db 



opt 
1 



s^P^V^a ASV^a 
It corresponds to numerical results (Table UTTll . 



(196) 
(197) 
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We use equations (l4Tll . (I48ll to calculate the gain.: 



We get : 



1 aV b , , 

It is in very good agreement with numerical data (Table IIIII) . Gain can be very large if 
density is low. 

B. Ballistic mode in one dimension : exact result (cf. section IlII C 2j) 



= ^^(71+72 + 73) (200) 

7. - '^^^ (201) 

72 = ]^3/2^^^ [94{-l)e--^ - g^il)) U(-l)e'^ + (73(l)e'^J (202) 

den = 5(i(l)e ^1^2 + 5(i(-l)e -^1^2 + 512(1)6 ^1^2 + g^[-i)e^ l^l^ (203) 

_ L2 numinum2 , . 

2a 

numi = /i(l) + /i(-l)e ^2 + + + ^3 + ^4 (205) 

/i(e) = 2 (a^74(e) (/^ + i^i) + (eLa - Li)) (206) 

(Ti = (/2(-l) + /4(1, 1) + /3(1)) (207) 

^2 = (/2(1) + /4(1, -1) + /3(-l)) (208) 

a3 = ihi-l, 1) + /5(1) + /6(1)) e'^T^ (209) 

(74 = (/4(-l, -1) + /5(-l) + /6(-l)) e"'^T^ (210) 

/2(e) = (v^ + eLa) L2 (Li - L2) (73(e) (211) 

/3(e) = -L^Liv^ (/i + Li) (V^ + eL2 + eL^) (212) 

/4(ei, 62) = ha{h + Li) ((2L2 + eiLi) (eav^ + L2) + L?) (213) 

/5(e) = -eh^L2 (Li + L2) (2 (ev^ + L2) L2 + L^) (214) 
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/6(e) = LlLi (L2 (ev^ + L2) (Li + L2) - (/i + Li) (v^ + eLa - eLi)) (215) 

2a 2a 

num2 = ^1 + ^2 - 9A{l)e-^ (/7(1) + /8(1)) - ^74(-l)e-^ (^(-l) + /8(-l)) (216) 

e^ + e^J (217) 

?2 = 2L2 (/I (/i + Li) a - f - 1 (218) 

/7(e) = (a + (Li + eLs) /i) f e' ^1^2 + e' ^1^2 j (219) 

/8(e) = [eha + + eL?) |^-e^ ^ + (220) 

rfem = (ei + 6 + ^3 + ^4 + ^5 + ^6 + ii) (221) 

{1 = 2Li (/i + Li)a (222) 

6 = I^2V^ ((« + sinh + 2L2v^cosh (^)) (223) 

^3 = ^1^2 [oi sinh (^^^ - 2L1L2 cosh ^ (224) 

^4 = -L2^fa [a + 2Li/i + L^) cosh sinh (^^^^ (225) 

^5 = -2 (Li (/i + Li) a + L^) cosh cosh (226) 

^6 = -Ua {2h + Li) sinh (^^^ cosh (^^^) (227) 

= - ((2/i + Li) a + Ll) sinh (^^^ sinh (^^) (228) 

ciera2 = ^i(l)e ^1^2 +^^(-1)6^^1^2+^2(1)6 ^1^2 + ^2(-l)e^i^2 (229) 

^?i(e) = L2(Li + L2) (v^-eL2) (230) 

^2(e) = (2/i + Ll) {-e^/^ - + eLi) + ea^'/^ + _ ^^3 ^231) 

^3(e) = /iV" (-ev^-Ls) +e2L2Li (232) 

g^{e) = ((eL2 - Li) /i + (233) 

h = b-a (234) 

a = + (235) 

Ll = vin (236) 
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L2 = Vt2 (237) 

This result have been checked by numerical simulations and by comparison with known 
limits. 



C. Static mode in 3 dimensions : more comparisons between the analytical ex- 
pressions and the simulations (cf. section [VAU) 
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Figure 23: (Color online) Static mode in 3 dimensions. Study of the minimum : its location in 
the ri,r2 space (a), and its value (b). sim means values obtained through numerical simulations, 
th means analytical values. Value expected if there was a prefect agreement between theory and 
simulations (black line), and values taking into account the simulations noise (dotted black lines) 
(we performed 10 000 walks for each point), a = 0.01 (brown squares), a = 0.1 (red crosses), a = 1 
(purple circles), a = 10 (blue stars), a = 100 (green diamonds). V = 1, k = 1. 

The numerical study of the minimum mean search time (Figl23l) shows that the analytical 
values gives the good position of the minimum in ti and t^, as soon as h/a is not too small. 
However, the value of the minimum is underestimated by about 10%. 

D. Diffusive mode in 3 dimensions 

1. Full analytical expression of tm (cf. section W B ^) 



1 



b^a^dpD 
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(X + Y + Z) 



(238) 



with 



(239) 

Y = 3^^ (240) 

Z=-2/30^ ^ '-^ ^ (241) 

(X 

a = ^/{tiD)-^ + {t^D^)-^ (242) 

D2 = ^V't2 (243) 

DD2 , , 

c^P-^ (244) 

"2 = (r2/^2)~' (245) 

a 6 tanh (a (6 — a) ) — 1 

^ = ^ \ \ n 246 

a — tann [a [h — a)) 

^ {a%a — 1) tanh {a {h — a)) + a {h — a) (247) 
ah — tanh (a (6 — a)) 

TT = ^ (248) 
tanh ( ^2 a) 



^. Dependence of tm with ti 

The mean detection time is very weakly dependent on ti as long as ti < 6D/v^ (Fig l24l) 



5. tm ^/*e regime of diffusion alone (cf. section VB_4_) 



We take a diffusive random walk starting from r = ro in a sphere with reflexive boundaries 
at r = 6 and absorbing boundaries at r = a, we get the following equation for t(ro) the mean 
time of absorption : 



1 / d 



Deffir, ( i- ( rr-^ ) ) = -1 (249) 



With the boundary conditions, the solution is : 



dt{ro)\\ ^ 
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Figure 24: (Color online) Diffusive mode in 3 dimension, tm from (I238p . tm(b ^> a,ri 0) from 
(fT44|l . T2 = r"^*'*'' (Il45|), D = 1,V = 1, a = 10 (dotted lines), a = 100 (lines), a = 1000 (symbols), 
b/a = 10 (blue, circles), b/a = 100 (red, squares). 



Than we average on ro, as the searcher can start from any point of the sphere v^ith the same 
probability : 



t 



1 



diff 



15Dab^ 



In the limit 6/a ^ 1 : 



b' 

3D^ 



(251) 



(252) 



4- Criterion for intermittence : additional figure (cf. section W B 5|) 



The figure [25] shovv^s the dependence of t^* v^^ith a. 



E. Ballistic mode in 3 dimensions 



1. without intermittence (cf. section IV C ^) 



In the regime without intermittence, ri is not necessarily 0. We calculate tm in two 
limits : Ti small or Ti large. 

a. Limit T2 —>■ 0, ViTi < a 
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Figure 25: (Color online) Diffusive mode in 3 dimension. Simulations : h/a = 2.5 (red crosses), 
b/a = 5 (blue squares), b/a = 10 (green circles), b/a = 20 (brown diamonds). Analytical expressions 
in the low target density approximation (6/a » 1) : ti = (|144p (with T2 = t2^^'^^ (|145p ) (dotted 
line), diffusion alone (|150p (continuous line). V = 1, D = 1 



In the limit viTi < a, we can consider phase 1 as diffusive, v^ith : 

D = ^tf n (253) 

We use the approached expression of obtained in the diffusive mode fll50p with this 
effective diffusive coefficient. 

= (5b^a^ + 56^ - 9b^a - a^) (254) 

5vfTia¥ ^ ' 

And in the limit b ^ a : 

_ b^ 

tm — 2 (255) 

ViTia 

b. limit r2 ^ 0, Ti ^ oo 

We name Voi the volume of the sphere. g{t) is the volume explored by the searcher after 
a time t. The volume explored during dt is nvio^dt. If we consider that the probability 
to encounter a unexplored space is uniform, which is wrong at short times but close to the 
reality at long times, the average of first explored volume at time t during dt is ^°'-~^^^*\ via?dt. 
Then in this hypothesis, g{t) is solution of : 

g[t) = / — Tivia du (256) 

Jo KjZ 
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This equation can be simplified taking a renormalized time r as r 



t, and / = g/Voi : 



;i - f{w))dw 



(257) 



Then, as /(O) = (nothing has been explored at time 0), /(r) = 1 — e~'^ . The probability 
to encounter the target at time t during dt (and not before) is the newly explored volume 
at time t divided by the whole volume Voi if we make the mean-field approximation. Then 
the probability p{r) than the target is not yet found at time r is solution of : 



dp 
dr 



(258) 



As p(0) = 1, the result is p{r) = e ^. Then the mean detection time of the target is 1 is 
renormalized time, is to say in real time : 

4^3 



'al 



?)a?vi 



(259) 



c. Numerical study 




0.05 0.1 0.5 1 5 10 50 



Figure 26: (Color online) Ballistic mode in 3 dimensions. Regime without intermittence {t2 = 0). 
\n.{tm/b^) as a function of In(ri), simulations for 6 = 5 (green squares) and 6 = 20 (blue circles). 
Ballistic limit (ri — > oo) (no intermittence) (|152p (red horizontal line), Diffusive limit {vti < a) 
(|254p with 6 = 5 (green dotted line), 6 = 20 (blue small dots), 6 » a limit (|255p (black line), a = 1, 
vi = 1. 



These expressions give a very good approximation of the values obtained through simu- 
lations (Figlini Fig l2T|) . In the regime without intermittence, tm is minimized for ri ^ oo. 
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Figure 27: (Color online) Ballistic mode in 3 dimension, vf as a function of ln(6), through simula- 
tions, a = 1,V = 1. 



2. Numerical vf ( cf. section \V C4\ ) 

In simulations Figl27l when b is small decreases, but stabilizes for larger b, which is 
coherent with the fact that this value is obtained through a development in b. The value of 
vi for large b is different (even if close) to the expected value. The main explanation of this 
discrepancy is that in the intermittence regime, the approached value of tm is about 20% 
away from the value obtained through simulations. 
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